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Abstract 

The girth of a graph is the length of its shortest cycle. We give 
an algorithm that computes in 0{nlog^ n) time and 0(n) space the 
(weighted) girth of an Ji- vertex planar digraph with arbitrary real edge 
weights. This is an improvement of a previous time bound of 0(n^/^), 
a bound which was only valid for non-negative edge-weights. Our 
algorithm can be modified to output a shortest cycle within the same 
time and space bounds if such a cycle exists. 



1 Introduction 

The girth of an unweighted graph is the length of a shortest cycle in the 
graph or oo if the graph is acyclic. It is a well-studied graph characteristic 
and has been shown to be related to numerous other properties of graphs, 
including vertex degree, diameter, connectivity, maximum genus, and vertex 
colourings [2],|1]. For instance, it is easy to see that for a graph containing a 
cycle, its girth is at most twice its diameter plus one. 

The problem of computing the girth of a graph has received some atten- 
tion. An 0{mn) time algorithm is known where m and n are the number of 
edges and vertices, respectively [Ej. Finding the length of a shortest cycle of 
even length can be done in 0{n^) time [TT] . 

In this paper, we focus on the problem of computing the girth of planar 
graphs. For this class of graphs, faster algorithms are known. A linear 
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time algorithm is presented in [3] but it only applies when the graph has 
bounded girth. For general planar graphs, Djidjev gave an 0(n^/^ log n) 
time algorithm. This can be improved to 0(n log n) time by applying the 
minimum cut algorithm of |T] to the dual graph. Recently, it was shown how 
to find the girth in O(nlogn) time [TO] . 

All these results for planar graphs assume that the graph is undirected. 
For planar digraphs, Weimann and Yuster [TU] gave an 0(n^/^) time algo- 
rithm and they asked whether a faster algorithm exists. 

The girth of a graph is defined for unweighted graphs but this definition 
immediately extends to the case where edges have real weight. Of the above 
algorithms for planar graphs, only the O(nlog^n) time algorithm in [T] and 
the 0{n^^'^) time algorithm in [TU] can handle weighted graphs and only if all 
edge weights are non-negative. The 0(?7,log n) time algorithm only applies 
when the graph is undirected. 

We consider the most general version of the problem for planar graphs. 
We show how to find the girth of an n-vertex planar digraph with arbi- 
trary real edge weights (non-negative as well as negative) in O(nlog^n) time 
and 0{n) space. In particular, we answer the question by Weimann and 
Yuster [10] of whether an o{n^^'^) time algorithm exists for computing the 
girth of a planar digraph. Our algorithm can output a shortest cycle within 
the same time and space bounds, assuming such a cycle exists. 

The organization of the paper is as follows. In Section [21 we introduce 
some basic definitions and notation. We give our algorithm for computing 
the girth of a planar digraph in Section [3] and in Section HI we show how 
to extend this algorithm to output a shortest cycle. Finally, we make some 
concluding remarks in Section [S] 

2 Definitions and Notation 

For a graph H, we let Vh and Eh denote its vertex and edge set, respectively. 
Let G = {V, E) be a digraph with real edge weights defined by weight function 
w : E ^ M.. For vertices u,v & V, we let dG{u,v) denote the length of a 
shortest path in G from u to v w.r.t. w (we omit w in the notation but this 
should not cause any confusion). If no path from u to v exists, we define 
dciu, v) = oo and if there is a path from m to f containing a negative-weight 
cycle, dciu.v) = — oo. 

The (weighted) girth of G is the length of a shortest cycle in G w.r.t. w. If 
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G is acyclic, we define its girth to be oo. If G contains a negative- weiglit cycle 
C, cycles of arbitrarily large negative weiglit can be obtained by traversing 
G sufficiently many times so in this case, we define the girth of G to be —oo. 

3 Computing the Girth 

In the following, let G = (V, E) be an n- vertex planar digraph with real edge 
weights defined by weight function w : E ^ M.. In this section, we show how 
to compute the girth of G in 0(n log^ n) time and 0{n) space. 

We may assume that G contains no negative-weight cycles since the al- 
gorithm in [S] can detect such cycles within our time and space bounds. If a 
negative- weight cycle is present, our algorithm outputs — oo as the girth of 
G. 

We will assume that G is triangulated with pairs of oppositely directed 
edges. If it is not, this can be achieved by adding edges of sufficiently high 
weight W so that finite shortest path distances in G will not be affected. We 
define 

where E' is the set of edges in the original graph. This way, we avoid dealing 
with infinite shortest path distances. And we can still detect the case where 
the girth is oo in the original graph since this holds if and only if the girth 
of the triangulated graph is at least 1 + YlieeE' 1^(^)1- 

We will identify G with a fixed plane embedding of the graph. 

Next, we apply the linear time cycle separator theorem of Miller [H] to 
G. This gives a simple cycle C in G (ignoring edge orientations) containing 
0{y/n) vertices such that at most 2n/3 vertices of G are in the closed bounded 
region Ri resp. closed unbounded region R2 of the plane defined by G . Let 
Gi resp. G2 be the subgraph of G containing the set of vertices and edges 
of G in i?i resp. i?2- If an edge of G belongs to both Ri and i?2, i-e., to C, 
then we only add it to one of the two subgraphs, say Gi. This ensures that 
Gi and G2 are edge-disjoint. 

We recursively compute the girth gi of Gi and the girth (72 of G2- Let g 
denote the length of a shortest cycle in G that contains at least two vertices 
of G. Any simple cycle in G that is neither fully contained in Gi nor in G2 
must contain at least two vertices of G. Thus, T[vin{gi, g2, g} is the girth of 
G so let us consider the problem of computing g. 
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Let H be the complete digraph with vertex set Vh = Vc and with weight 
function wh ■ — > M defined by wh{u,v) = mm{dGi{u,v),dc2{u,v)} for 
all distinct u, v G Vh- 

For any u,v & Vh, a shortest path from n to f in G can be decomposed 
into subpaths each of which has the property that it is a shortest path in 
either Gi or in G2 with both its endpoints on C. It follows that dniu^v) = 
dciu^v). 

We apply the algorithm in [7] to compute da^iu^v) and dciiu^v) for all 
G Vc using a total of O(nlog^n) time over all recursion levels. Hence, 
we obtain H and its edge weights in this amount of time over all recursion 
levels. 

Any shortest cycle in G containing at least two vertices of G can be 
decomposed into subpaths each having the property above. Hence, such a 
cycle has the same length as a shortest cycle in H so H has girth g. 

What remains therefore is to find the length of a shortest cycle in H. We 
will show how to do this in 0{n\og^ n) time. Since there are 0(log?T,) recur- 
sion levels, this will imply that the girth of G can be obtained in 0{n\og^ n) 
time. 

Let Ml, . . . , Mm be the vertices of G. To compute the girth g of H, we will 
compute, for i = 1, . . . , m, the length gi of a shortest cycle in H containing 
Mj. Then it is clear that we can obtain the girth of H a.s g = minj^fi, . . . , gm}- 

We first reduce our problem to one where all edge weights are non- 
negative. This is done as follows. We compute single source shortest path 
distances with source, say, mi, using the 0(?T,log^n) time Bellman- Ford vari- 
ant of [5] (in fact, it can be done in only 0{na{n)) time with ideas from [8] 
but this will not improve our bound). 

Then in 0[\Eh\) = 0{n) additional time, we can obtain the reduced cost 
w~l^{e) of each edge e = (ui,Uj) of H (w.r.t. wh)'- 

w^{e) = dniui, Ui) + Wnie) - dniui, uj). 

By the triangle inequality, > and it is easy to see that for any cycle in 
H, its length w.r.t. wh is identical to its length w.r.t. w~^. This gives us the 
desired reduction. 

For each pair of vertices Ui and Uj in H, we let dfj{ui,Uj) denote the 
shortest path distance from Ui to Uj in H w.r.t. w^. 

Now, let us consider the problem of computing gi for some i. Since 
gi = mm{d'jj{ui,Uj) + wjj{uj,Ui)\j = l,...,m,j 7^ i}, we can obtain this 
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value in 0{m) = 0{-/n) time in addition to the time for computing single 
source shortest path distances from Ui w.r.t. w~^. And since > 0, we can 
apply the Dijkstra variant of [5] to compute these distances in 0(m log^ n) = 
0(-\/nlog^n) time. Over all i, this is O(nlog^n), as desired. 

It follows from the above that we can find the girth of H in O(nlog^n) 
time and we can conclude this section with the following result. 

Theorem 1. The girth of an n-vertex planar digraph with real edge weights 
can be computed in O(nlog^n) time and 0{n) space. 

Proof. We have shown the time bound above and since the algorithms in [5l 
[71 [HI [9] all require 0{n) space, this bound also holds for our algorithm. □ 

4 Finding a Shortest Cycle 

We now show how to extend the algorithm of the previous section to compute 
a shortest cycle in G within the same time and space bounds. We may assume 
that G contains no negative-weight cycles since otherwise, a shortest cycle 
does not exist. 

With the definitions above, it suffices to find in 0{n log^ n) time and 0{n) 
space a shortest cycle in G containing at least two vertices of C . As already 
noted, such a cycle corresponds to a shortest cycle in H . The algorithm 
above finds vertices Ui and Uj of C such that a shortest path Pij in H from 
Ui to Uj followed by the edge {uj,Ui) is a shortest cycle C' in H. The edges 
on Pij can be found by traversing the shortest path tree in H rooted at Ui. 
Since this shortest path tree has already been computed, we can thus find 
all edges of C' within the required time and space bounds. 

What remains is to replace each edge e of C' by a simple shortest path 
Pe in G between its endpoints. We will show how to do this in O(nlogn) 
additional time. 

To obtain these paths, we need to take a closer look at the multiple-source 
shortest path algorithm of Klein [7] which we applied to find the weights of 
edges of H. For Gj, i = 1,2, his algorithm maintains a dynamic tree data 
structure, which is initally a shortest path tree in Gi rooted at, say, Ui, then 
at U2, and so on. To obtain the weights of edges of H, this data structure is 
repeatedly queried, first for the distance in Gi from ui to all other vertices 
of C, then from U2, et cetera. 
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Now, we also need the actual paths in the shortest path trees correspond- 
ing to these distances. Once the edges of C have been identified, we can 
run Klein's algorithm again. During the course of this second run of the 
algorithm, we can query the dynamic tree data structure to obtain the paths 
corresponding to edges of C . Querying for a single vertex on a path takes 
O(logn) time so the total time is 0(A;logn), where k is the total number of 
vertices (with repetitions) over all paths. 

For this strategy to work, k should not be too big. We will ensure this by 
modifying Dijkstra without increasing its time and space bounds such that 
C will have the least number of edges among all shortest cycles in H. We 
then show that none of the paths defining edges of C will share any edges 
of G, implying that k = 0{n) and hence that a shortest cycle in G can be 
found in 0(nlog^n + klogn) = O(nlog^n) time and 0{n) space. 

To output a shortest cycle in H with the minimum number of edges, we 
modify the Dijkstra algorithm in [5]. We will not go through all the details 
but assume that the reader is familiar with the paper. 

During the course of that algorithm, we maintain not only shortest path 
distances w.r.t. wjj- but also w.r.t. the unit weight function. Now, whenever 
there is a tie between which vertex to pick next from the heap, we pick one 
for which the number of edges on a shortest path from the root of the par- 
tially built shortest path tree to Ui is minimized. This can be incorporated 
into the algorithm without increasing its time and space bounds by lexico- 
graphically ordering vertices, first according to weighted and then according 
to unweighted distances from the root of the shortest path tree. 

It remains to show that if C is a shortest cycle in H and it is picked such 
that it has the minimum number of edges then none of the shortest paths in 
G corresponding to edges of G' share edges. 

So let Pi and P2 be shortest paths in G corresponding to distinct edges 
ei = {ui,Uj) and 62 = {ui',Uj') in H, respectively. Assume for the sake of 
contradiction that e = {u,v) is an edge shared by Pi and P2. Let P be the 
subpath of G' starting in uj and ending in Ui'. Without increasing the length 
of G', the subpath P1PP2 of G' can be replaced by the path P' defined as 
the prefix of Pi ending in v followed by the suffix of P2 starting in v. Since 
Gi and G2 are edge- disjoint, either Pi and P2 both belong to Gi or both 
belong to G2- Hence, P' is a path in either Gi or in G2 so we can replace ei 
and 62 by the edge (uj, Uj>) in H. But this reduces the number of edges of 
G' without increasing its length, contradicting the choice of the cycle. 

It follows that none of the shortest paths in G corresponding to edges of 
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C share edges. By the above, this suffices to show the following. 

Theorem 2. A shortest cycle in an n-vertex planar digraph with real edge 
weights can be computed in 0{n\og^ n) time and 0{n) space, assuming such 
a cycle exists. 

If G has girth — oo, a shortest cycle does not exist. But we can still output 
a negative-weight cycle within our time and space bounds by applying the 
algorithm in [S]. 

5 Concluding Remarks 

We showed how to compute the girth of an ra-vertex planar digraph with 
real edge weights in O(nlog^n) time and 0{n) space. This is a significant 
improvement over the previous best bound of 0(n^/^) which only applied 
to planar digraphs with non-negative edge weights. We also showed how to 
output a shortest cycle without an increase in time or space, assuming such 
a cycle exists. 

In |5], it is suggested that the results of that paper can be generalized 
to the class of bounded genus graphs. If this is true, we believe that a 
generalization of our algorithm to this class is also achievable. 
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